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Strong Disorder Renormalization is an energy-based renormalization that leads to a complicated 
renormalized topology for the surviving clusters as soon as d > 1. In this paper, we propose to 
include Strong Disorder Renormalization ideas within the more traditional fixed cell-size real space 
RG framework. We first consider the one-dimensional chain as a test for this fixed cell-size procedure: 
we find that all exactly known critical exponents are reproduced correctly, except for the magnetic 
exponent (3 (because it is related to more subtle persistence properties of the full RG flow). We 
then apply numerically this fixed cell-size procedure to two types of renormalizable fractal lattices 
(i) the Sierpinski gasket of fractal dimension D = In 3/ In 2, where there is no underlying classical 
ferromagnetic transition, so that the RG flow in the ordered phase is similar to what happens in 
d — 1 (ii) a hierarchical diamond lattice of fractal dimension D = 4/3, where there is an underlying 
classical ferromagnetic transition, so that the RG flow in the ordered phase is similar to what 
happens on hypercubic lattices of dimension d > 1. In both cases, we find that the transition is 
governed by an Infinite Disorder Fixed Point : besides the measure of the activated exponent tp, we 
analyze the RG flow of various observables in the disordered and ordered phases, in order to extract 
the 'typical' correlation length exponents of these two phases which are different from the finite-size 
correlation length exponent. 
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INTRODUCTION 



The choice to work in real-space to define renormalization procedures, which already presents a great interest for 
pure systems [![, becomes the unique choice for disordered systems if one wishes to describe spatial heterogeneities. 
. Whenever these disorder heterogeneities play a dominant role over thermal or quantum fluctuations, the most ap- 
propriate renormalization procedures are Strong Disorder renormalizations that have been introduced by Ma and 
Dasgupta Q : as shown by Daniel Fisher [1, these Strong Disorder RG rules lead to asymptotic exact results if 
the broadness of the disorder distribution grows indefinitely at large scales. In dimension d = 1, exact results for a 
large number of observables have been explicitly computed 0, [B| because the renormalized lattice of surviving degrees 
■ of freedom remains one-dimensional. In dimension d > 1, the Strong Disorder RG procedure can still be defined, 
J> ' but it cannot be solved analytically, because the topology of the lattice changes upon renormalization. Nevertheless, 
. Strong Disorder RG rules can be implemented numerically. For instance, for the disordered Quantum Ising model, 
these numerical RG studies have concluded that the transition is also governed by an Infinite Disorder fixed point in 
dimensions d — 2, 3, 4 in agreement with the results of independent quantum Monte-Carlo in d = 2 [l7l Il8j|. 

Nevertheless, the complicated topology that emerges between renormalized degrees of freedom in dimension d > 1 
tends to obscure the physics, because a large number of very weak bonds are generated during the RG, that will 
eventually not be important for the forthcoming RG steps. In this article, we propose to include strong disorder 
RG ideas within the more traditional fixcd-length-scale real space RG framework that preserves the topology upon 
renormalization. In particular for renormalizable fractal lattices like the Sierpinski gasket or diamond hierarchical 
lattices, this fixed-length-scale RG procedure allows to use the so-called 'pool method' and to study numerically very 
large system sizes. (Note that when the full strong disorder renormalization is applied to fractal lattices like the 
Sierpinski gasket as in Ref [l9j . the system sizes are limited because one cannot take advantage of the self-similarity 
of the original lattice which is immediately broken by the full RG.) 

The paper is organized as follows. In section UH we introduce the principles of Strong Disorder RG within a fixed 
cell-size RG framework. In section [TTT1 we present our numerical result for d — 1, and conclude that this procedure 
correctly captures all critical exponents except for the magnetic exponent /? which is related to persistence properties 
of the full RG flow. In section IIV[ we apply numerically our procedure to the Sierpinski gasket of fractal dimension 
D = In 3/ In 2 = 1.58.. with no underlying classical ferromagnetic transition. In section |Vj we study numerically 
the case of a diamond hierarchical lattice presenting an underlying classical ferromagnetic transition. Section IVII 
summarizes our conclusions. Appendix lAl contains a reminder on the usual Strong Disorder RG on arbitrary lattices 
and on the properties of renormalized observables as a function of the energy-based RG scale V. 
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II. STRONG DISORDER RG WITHIN A FIXED CELL-SIZE REAL-SPACE RG FRAMEWORK 

In this paper, we consider the quantum Ising model defined in terms of Pauli matrices 

« = - E JijViVj -J2 h ^ w 

<i,j> i 

where the nearest-neighbor couplings Jjy > and the transverse fields hi > are independent random variables 
drawn with two distributions ^coupling (J) an( i ^fieidQ 1 ). 

As recalled in Appendix [XJ the Strong Disorder Renormalization for the quantum Ising model of Eq. [1] is an 
energy-based RG , where the strongest ferromagnetic bond or the strongest transverse field is iteratively eliminated. 
In dimension d > 1 this leads to a complicated topology for the network of surviving clusters at large RG scale T. In 
this section, we propose to apply the Strong Disorder RG principles within a fixed cell-size framework , i.e. the RG 
scale will not be an energy-based scale like T, but a length-scale L as in usual real-space RG procedures. Let us first 
describe the expected scalings as a function of L in various phases, to justify the possibility of such a fixed-length 
procedure. 

A. Expected scaling of renormalized observables as a function of the size L 

Since the energy-RG-scale T is associated to some length-scale It via Eq. IA91 it seems clear that all RG flows as a 
function of T described in Appendix [S] can be reformulated as RG flows as a function of the length scale L. Physically, 
these scalings in L simply describe the scaling of observables of finite-size samples. We should stress however that 
the change from the energy-scale T to the length-scale L is not just a 'mathematical change of variables', because 
probability distributions are involved. Indeed, the ensemble of disordered samples of volume L d is characterized 
by a probability distribution PL^iast) of the last RG scale Ti ast - Similarly for a sample in the thermodynamic 
limit L — > +oo in d = 1, the state at RG scale V is characterized by probability distribution of lengths (Ib,Ic) for 
renormalized bonds and renormalized clusters Q (in dimension d > 1, the state at RG scale T involves in addition a 
complicated topology of the renormalized surviving clusters) . 

1. Critical Point : RG flow as a function of the size L 

At the Infinite Disorder critical point described in section IA 41 the renormalized transverse fields and the renor- 
malized couplings remain in competition at all scales, i.e. they display the same activated scaling in L with some 
exponent tjj 

lnJ L =-L v V 

kih L =-Vi'v c (2) 

where {u c ,v c ) are 0(1) random variables. The exponent ip corresponds to the activated exponent of the gap G(L) oc 
e- L * fEa lATl . 

2. Critical Region : Finite-Size Scaling with the exponent vfs 
At this Infinite Disorder fixed point, one expects the following finite-size scaling form for the typical values 0] 

In f L vp = huh = -L'^Fj (l 1 '^^ - h c \) 

\xih l l p = WhZ = -L^F h (L x ^ ps \h - h c \) (3) 
as well as for the width of the distribution of (In Jl) 

1 /2 

A L = ((hTjIF - (InX) 2 ) = L*F A (h 1 '"^ \h ~ h c \) (4) 

(and similarly for the width of distribution of (ln/i^)). Here vp$ is introduced as the correlation length exponent 
that govern all finite-size effects in the critical region. 
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The magnetization of critical clusters grows with the fractal dimension df introduced in Eq IA21I The intensive 
magnetization is expected to follow the usual power-law finite-size scaling form 

m = cx L~ x Fm [l^ Ufs \h - h c \j with x = d-d f (5) 
3. Disordered Phase : RG flow as a function of the size L 



In the disordered phase described in section [A 31 the transverse helds hi are not renormalized anymore asymptoti- 
cally, i.e. they converge towards finite values as L — > +oo 



\nh t yP = \nh OD (x-(h-h c )- K (6) 
where the exponent k of the essential singularity satisfies 

K = tjjVFS (J) 

as a consequence of the matching with the finite-size scaling form of Eq. [31 

The renormalized couplings Jl is expected to have the same scaling as the two-point correlation function 0, 0] 

In J L = -— + L UJ A(h)u (8) 

The first term is non-random and describes the exponential decay with the size L, where £t yp represents the typical 
correlation length 

\nj L yp = huE oc -J- (9) 

The compatibility with the finite-size scaling form of Eq. [3J implies that the typical correlation exponent u typ is 
different from the finite-size correlation length exponent vps and reads 

v tyv = (1 - iP)vfs (10) 

The exponent vps is expected 0, 0] to correspond to the exponent v av of the averaged two-point correlation function 
(v FS =u av ). 

The second term in Eq[8]contains an 0(1) random variable u. It is usually subleading with respect to the first term, 
i.e. of order with some exponent u> < 1. We have argued in [2(| that this exponent to should coincide with the 
droplet exponent u>dp(D = d — 1) of the Directed Polymer with D = (d— 1) transverse directions. The compatibility 
with the finite-size scaling form of Eq. [4] for the width Al of the distribution of In Jl 

1/2 



A L = f(lnJ L ) 2 -(lnJ L ) 2 ) cx A{h)L" (11) 

\ / L— ¥00 



implies that the amplitude A(h) is non-singular as h h c only if ip = u>, which is known to be the case in d = 1 [j. 
On the other hand, if the two exponents turn out to be different ip ^ oj, then the amplitude A(h) has to present the 
following power-law singularity 

A(h) oc (h-h c )-W~ u)v * s (12) 

h^h+ 

4- Ordered Phase : RG flow as a function of the size L 



In the ordered phase described in section IA 51 the magnetization /i of surviving clusters grows extensively 

fX L cx L d (13) 

i.e. the intensive magnetization is finite and vanishes with some power-law singularity compatible with the finite-size 
scaling form of Eq. [5] 

m = ^4 cx (h c — h)P with (3 = xups (14) 
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Since surviving clusters have an extensive magnetization, the logarithm of their renormalized transverse fields grows 
also extensively 

ln/iL oc - (^) v (15) 

X->oo \£h/ 

where v is an 0(1) random variable. The length scale £^ represents the characteristic size of finite disordered clusters 
within this ordered phase. The compatibility with the finite-size scaling form of Eq. |3] implies that the typical 
correlation exponent Vh reads 

Vh = (l ~ ^) vfs (16) 

Note that plays in the ordered phase a role similar to v tyv in the disordered phase, but that they coincide only if 
d= 1 (EqEEDJ). 

For the asymptotic behavior of the renormalized couplings Jl between surviving clusters, one needs first to determine 
whether there exists an underlying classical ferromagnetic transition, as recalled in section [A 51 

(a) When there is no underlying classical ferromagnetic transition, the typical renormalized coupling remains finite, 
and presents the same essential singularity as in Eq. [6] 



\nJZ v = \^Joo <x-(h c -h)~ K (17) 

(b) When there exists an underlying classical ferromagnetic transition, the renormalized couplings Jl do not remain 
finite as in Eq. [T7] but grow asymptotically at large L with the scaling of the classical random ferromagnet model 

J L oc (TL d °+L e °u (18) 

L— >+oo 

The first non-random term grows as the surface L ds of dimension d s (a being the surface tension) . The second term 
contains a random variable u and grows with some exponent 8 S . For instance in d = 2 where the surface of dimension 
d s = d — 1 = 1 is actually a line, the exponent 9 S coincides with the droplet exponent ojdp — 1/3 of the Directed 
Polymer in dimension 1 + 1. In dimension d = 3 where the surface has dimension d s = d — 1 = 2, the exponent O s 
has been numerically measured to be Q s ~ 0.84 [21]. It is clear that in the regime of large L where Eq. [TBI holds, 
the strong disorder RG procedure is not appropriate anymore, because the couplings J grow with the scale L (so 
the decimation of the biggest coupling tend to create a bigger renormalized coupling) , and because the width of the 
probability distribution of (In Jl) actually decrease with L (instead of becoming broader and broader) 



A, 



f(lnJ L )2-(lnJ L ) 2 ) 1/2 oc L- d '+°' (19) 



This is why in the language of the energy scale T, the strong disorder RG stops at some finite value T perco where 
percolation occurs. In the language of the length scale L, the RG flow exists for all L but is non-monotonous in the 
critical region of the ordered phase : at the beginning, the typical coupling decays in Eq. [2] and Ax, grows as lr , 
whereas asymptotically the typical coupling grows as in Eq. [18] and A l decays as in Eq. [T9l 

5. Discussion 

The fact that RG flows can be reformulated in terms of the length-scale L suggests that some fixed cell-size real- 
space RG procedure should be possible even for Infinite Disorder fixed points. In the following, we propose such an 
explicit RG procedure, where Strong Disorder decimation rules are used within a fixed cell-size framework. 

B. Principles of the fixed cell-size RG procedure for the one-dimensional chain 



For clarity, we first explain the ideas on the case of the one-dimensional chain, and compare with the exact results 
of the energy-scale strong disorder RG 
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1. Notations 



Let us first introduce notations. At generation n, corresponding to a length L n = b n , a renormalized open bond 
]A, B[ is characterized by the correlated renormalized variables ( Jab, r A, Aha', tb, Ahb) distributed with some joint 
probability distribution P n ( Jab, t a, Aha', r B, Ahb) '■ 

(i) Jab represents the renormalized ferromagnetic coupling between the two end points A and B 

(ii) rA represents the multiplicative factor that comes from the interior of the bond and that should be applied to 
the transverse field of A, see the RG rule of Eq. A3. (Similarly tb concerns the other end B of the bond). 

(iii) A ha represents the excess of magnetization that comes from the interior of the bond and that should be added 
to the magnetization of A, see the RG rule of Eq. A4. (Similarly Ahb concerns the other end B of the bond). 

At generation n = : the coupling Jab is just an initial coupling of the model drawn with some disorder distribution 
^coupling (Jab); the sites A and B have magnetic moments ha = 1 = M-B, corresponding to Aha = = Ahb ', and 
the transverse fields and hjp are the initial transverse fields of the model, corresponding to ta = 1 = r B . So the 
initial joint-distribution is simply 

Pu=o(Jab , r A, Aha', r B, Ahb) = n coup ii n g(J AB )8(rA - l)8{AH A )8{r B - I)S(Ahb) (20) 



2. Step 1 

The first step consists in building the elementary structure of generation (n+l) from the renormalized open bonds 
of generation n. 

For the one-dimensional chain with a scaling factor b, the explicit construction is as follows : 

la) Draw b independent renormalized open bonds at generation n, where each open bond j = 1, 2, .., b is characterized 
by the variables (Ja 6 b 3 , ta 6 , AhAj ; r B] , Ahb 3 ■ 

lb) For each j — 1,2, ..(b — 1), connect the two open bonds (j) and (j + 1) in series via the introduction of 
an intermediate site Cj corresponding to the merging of Bj with Aj + i. The intermediate site Cj has thus for 

l of its own initial magnetization He] = 1 a 
the renormalization of the interiors of the two neighboring bonds 



magnetization the sum of its own initial magnetization He — 1 an d of the excesses Aha j+1 and Ahb^ coming from 



Mc 3 = + ^a ]+1 + Ahb, (21) 

Similarly its transverse field is the product of its own initial transverse field h^} and of the two possible reducing 
factors rA j+1 and coming from the renormalization of the interiors of the two neighboring bonds 

h C] = fcg r Aj+1 r B] (22) 

The end-point A of the new bond of the (n + l) generation will correspond to the end-point A\ of the bond j = 1, 
and inherits the associated variables 

AH { f epl) =AH Al (23) 

and 

r { f epl) =r Al (24) 

Similarly, the end-point B of the new bond of the (n + l) generation will correspond to the end-point Bb of the bond 
j = b and inherits the associated variables 

Anf ep 1} - Ahb, (25) 



and 

r Bb (26) 



.(step 1) 

' 'rib 
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3. Step 2 

The second step consists in applying the usual strong disorder RG rules (recalled in section IA 1J to the internal 
structure with renormalized bonds of generation n that has been constructed in step 1, up to the final state containing 
a single renormalized open bond of generation (n + 1). 

For the one-dimensional chain with a scaling factor 6, the explicit strong disorder renormalization of the internal 
structure is as follows : 

These two end-points A and B of the new open bond of the (n + 1) generation are non-decimable by defini- 
tion, whereas all (6— 1) intermediate points (C\, C%, .., C&_i), and all internal bonds (Jaci — Ja 1 b 1 ,Jc 1 c 2 — 
Ja 2 b 2 , Jc b - X B = JA b B b are decimable. 

We apply to this internal structure the usual strong disorder RG rules recalled in Appendix Al. So we look for the 
maximum between the 6 couplings and the (6—1) transverse fields 

Q = max [Jac x , hc ± , Jc t c 2 , hc 2 , hc h _i, Jc^b] (27) 

and we apply iteratively the strong disorder RG rules up to the full decimation of this internal structure : the final 
output is the the joint variables (Jab, r A, A[ia', t b, A/z#) of an open bond of generation (n + 1). 



4- Explicit RG for the rescaling factor b — 2 



Step 1: 



Step 2: 



J 



A l B l 



5, 



C (h c ,\l c ) 
a =msLx(J Ai B l ,J Al B 1 ,h c ) 



J. 



B 2 



A=A\ 



B= B ? 



'AB 



(r B ,A\i B ) 



FIG. 1: To construct an open bond of generation (n + 1) with its renormalized variables ( Jab, ta, A/^a; rn, A/xs) defined in 
section Hi B II there are two steps : 

Step 1 : one draws 6 (here b — 2) independent open bonds of generation n with their renormalized variables, and one connects 

them via intermediate points Cj with renormalized variables (see section Hi B 2 1 for more details). 

Step 2 : one applies the usual strong disorder RG rules to the internal structure (see section Hi B 31 for more details). 

The corresponding RG equation for b = 2 is written explicitely in Eq. [39] 

As an example, we now describe more explicitly the case of rescaling factor 6 = 2 (see Fig. [T]). Step 1 con- 
sists in drawing 6 = 2 independent renormalized open bonds at generation n, with their renormalized variables 
(J Al B 1 ,rA 1 , AnA^rsj,, A/UbJ and (J A2 B 2 ,rA 2 , AfJ-A 2 ;r B2 , A^s 2 ), and in the merging of A 2 and B x into a single 
internal site C that has thus for renormalized transverse field 



h c = n§' r A2 r Bl 



(28) 



and for renormalized magnetic moment 



V-c = He + + A/x Bl 



(29) 
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Step 2 consists in the Strong Disorder Renormalization of the internal structure made of the site C and of the two 
ferromagnetic bonds surrounding C, i.e. Eq. 1271 reduces to 



fi = max 



J AC = JAtB^hc = h i § ) r A 2 r B 11 JcB = Ja 2 b 2 (30) 



( i) If = he, then the site C is decimated according to the usual Strong Disorder RG rules of Eq A2 : the 
renormalized coupling reads 

Jab = (31) 
h c 

whereas the properties of the two ends are unchanged, i.e. the supplementary factors coming from this second step 
are 

(step 2) _ 1 
r A — 1 

A^ tep 2) = (32) 

(ii) If O = J ac, then the site C is merged with the end A according to the usual Strong Disorder RG rules of Eqs 
A3 and A4 : the properties of the end A get renormalized according to 

{step 2) _ he 
r A 



J AC 

ep 2) _ 

A 



A^ tep 2) = Mc (33) 



whereas the renormalized coupling is (Eq A5) 



Jab = Jcb (34) 



(iii) If ft = Jcb, then the site C is merged with the end B according to the usual Strong Disorder RG rules of Eqs 
A3 and A4: the properties of the end B get renormalized according to 



(step 2) _ flC 
B 



B j 

Jc 



whereas the renormalized coupling is (Eq A5) 



Aa4 S ' CP 2) = (35) 



Jab = J ac (36) 



The renormalized open bond of generation (n + 1) is thus characterized by these correlated final variables 
( J^g al , r A mal , A[i A mal ; r^ nal , A^ iTlai ), where the coupling J A g al has been computed either with (i),(ii),(iii), and 
where the variables associated to the end-points A and B contain the contributions of the two steps described above, 
i.e. the excess of magnetizations of the two boundaries read 

A ^nal = ^(step 1) + ^(step 2) 

A/4 W = A/4 s * ep 1] + A/4f ep 2) (37) 



while their multiplicative factors read 



final (step 1) (step 2) 

A — r A r A 



final _ (step 1) (step 2) /„o^ 
B " ' B ' B \°°) 

Equivalently, the joint distribution P„( J, ta, A^a, tb, A/j,a) evolves according to the RG equation 
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P n +i ( Jab , r A , A(j, A , r B , A/i B ) = J dJ Al Bl dr Al dAfi Al dr Bl dAfi Bl P n (JA lBl ,r Al , A^i Al , r Bl , A/i Bl ) 
J dJ A2 s 2 dr A2 dAfi A2 dr B2 dA[i B2 P n [Ja 2 b 2 , r a 2 , A^a 2 , r B . 2 ,Afi B . 2 ) 



dh c dfi c S{h c - hr Bl r A2 )S(nc - (1 + A^ Sl + A^ 2 ) 



5(Jab - J AiB\J A2B2) ^fi^ r ^ _ r J 4 1 )<5(A/i J 4 - A(j,A 1 )S(rB - r B2 )5{A^i B - A[i B ) 
h c 

{Ja iBi > 

5{J A b - Ja 2 b 2 )S(t a - r Al -^—)S(Afi A - (AjUj + fic))S{r B - r B2 )S(Afi B - A^ B2 ) 

JA1B1 



[6{h c > J AlBl )0(h c > J MB2 ) 

xi j Ja 1 b 1 Ja 2 b 2 ) x xf 

o{Jab 7 )o{r A - i 

he 

+9{J AlBl > h c )0(J AlBl > J A2 e 2 ) 

{J AB - J A . 2B2 )S(r A - r Al ° )S 

Ja iBi 

+0{Ja 2 b 2 > h c )9(J A . 2B2 > Ja 1 b 1 ) 

he 

S{Jab - JA 1 B 1 )S(r A - r Al )5(A^i A ~ A^ Al )5{r B - r B2 - )5(A/i B - (Afj, B2 + n c )) (39) 

J,4 2 B 2 

where the two first lines corresponds to the draw of two independent open bonds of generation n, the third line 
corresponds to the merging of A2 and B\ into the intermediate site C, and where the last lines correspond to the 
three possible decimations (respectively of he, of J AlBl — J A c and of Ja 2 b 2 = Job)- 

C. Discussion 

We have described in detail the case d = 1. It is clear that in the limit of infinite block size b —> +00, one recovers 
the full usual Strong Disorder RG. In the next section Hill we have thus chosen to study numerically the 'worst' case 
of rescaling factor b = 2 to compare with the exact results of the full Strong Disorder RG corresponding to b — >• +00. 
The results with other block sizes b = 3, 4, ... are expected to interpolate between b — 2 and b —> +00. 

More generally, besides this one-dimensional case, we can apply the same strategy to other exactly renormalizable 
lattices. As discussed in section IA 5[ the properties of the ordered phase depend on the existence of an underlying 
classical ferromagnetic phase. We have thus chosen to consider below the two following cases : 

(1) in section HVl we consider the Sierpinski gasket where there is no underlying classical ferromagnetic phase. 

(2) in section [Vj we consider the hierarchical diamond Lattice where there is an underlying classical ferromagnetic 
phase. 

III. NUMERICAL RESULTS IN d = 1 FOR THE RESCALING FACTOR 6 = 2 

We have applied the procedure explained in detail in the previous section for the rescaling factor 6 = 2. In this 
section, we describe our numerical results and compare with the exact solution of the full Strong Disorder RG [H . 

A. Numerical details 

We have used the fiat initial distribution of couplings on the interval < J < 1 

^coupling 

(J) = 0(0 < J < 1) (40) 

and a uniform transverse field h. This choice of box-distribution is very common in numerical studies of random 
transverse field Ising models In addition in dimension d = 1, it is exactly known that any infinitesimal disorder 

flows towards the same Infinite-Disorder fixed point Q. 

The corresponding exact transition point given by the criterion In h = In J [4( is 

h exact = e W7 = e -i = 0.36787944... (41) 

We have performed numerical simulations with the so-called 'pool-method' which is very much used for disordered 
systems on hierarchical lattices, in particular for spin-glasses (22|, but also for disordered polymers [23|,[24j]. The idea 
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is to represent the joint probability distribution P n (J, Ta, A/^, tb, A/i^) at generation n, by a pool of M realizations 



.0 



A/i^, r|' , A/i^) where i = 1,2,...,M. The pool at generation (n + 1) is then obtained as follows : each 
new set ( J0, , , , A/x^) of generation (n + 1) is obtained by choosing & = 2 values i at random with their 
corresponding variables in the pool of generation n. 

The results presented in this Section have been obtained with a pool of size M = 10 7 . We find that the corresponding 
transition point satisfies 



.0 



0.367910 < h p c ool (M = 10 7 ) < 0.367912 



(42) 



i.e. it is slightly larger than the exact value of Eq. @J] as a consequence of the finite size M of the pool. We now 
describe our numerical results concerning the RG flow of the renormalized variables in various phases using < n < 50 
generations corresponding to lengths 



1< L = 2 n < 2 



50 



(43) 



B. RG flow of the typical renormalized coupling J 



' yp 




20 



In A, 



10 













*2?"h=0.367912 




h=0.39 


h=0.367910 


(b) 










h=0.32 





In L 



10 



20 



30 



In L 



40 



FIG. 2: (d = 1 ) (a) RG flow of the logarithm of the typical renormalized coupling (In j£ yp ) in a log-log plot : the slope is in 
the ordered phase (e.g. h — 0.32), the slope is 1 in the disordered phase (e.g. h = 0.39), and the slope is ^ — 0.5 at criticality 
(before the bifurcation of the curves corresponding to h — 0.367910 and h — 0.367912 ). (b) RG flow of the width Al of the 
distribution of the logarithm of the renormalized couplings in a log-log plot : the slope is in the ordered phase (e.g. h — 0.32), 
the slope is ui ~ 0.5 in the disordered phase (e.g. h — 0.39), and the slope is tp ~ 0.5 at criticality (before the bifurcation of the 
curves corresponding to h = 0.367910 and h = 0.367912 ). 



On Fig. [2] (a), we show our data concerning the RG flow of the typical coupling 



J 



typ _ „ln J L 



(44) 



Our results for the L-dependence in the two phases and at criticality 



oc 

L— >+oo 



Cst 

—L^ with V' 
-L 



0.5 



(45) 



are in agreement with the exact solution 



10 



C. RG flow of the width Al of the logarithms of the renormalized couplings 



On Fig. [2] (b), we show our data concerning the RG flow of the width of the distribution of the logarithms of 
the couplings. Again, our results for the L-dependence in the two phases and at criticality 



0.5 

0.5 (46) 



are in agreement with the exact solution 
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FIG. 3: ( d = 1 ) (a) RG flow of the logarithm of the typical renormalized transverse field (ln/i£ ,p ) in a log-log plot : the 
slope is d = 1 in the ordered phase (e.g. h — 0.32), the slope is in the disordered phase (e.g. h = 0.39), and the slope is 
ip ~ 0.5 at criticality (before the bifurcation of the curves corresponding to h = 0.367910 and h — 0.367912 ). (b) RG flow of 
the renormalized magnetization in a log-log plot : the slope is d — 1 in the ordered phase (e.g. h — 0.32), the slope is in 
the disordered phase (e.g. h — 0.39), and the slope is dj ~ 0.85 at criticality (before the bifurcation of the curves corresponding 
to h = 0.367910 and h = 0.367912). 



D. RG flow of the typical renormalized transverse field h? p 

On Fig. |3] (a) , we show the RG flow of the typical renormalized transverse field 



h tw = ( MhL ( 47) 

Again, our results for the L-dependence in the two phases and at criticality 

-L 

with V ^ 0.5 

(48) 

u— >™t-oo 

are in agreement with the exact solution 

E. RG flow of the renormalized magnetization \il 



In h^ p \ h<hc 


L- 


cx 

-> + <X) 




L- 


CX 

->+<x> 


lnh^ p \ h>hc 


L- 


cx 



On Fig. [3] (b), we show the RG flow of the magnetization /i_L of renormalized clusters as a function of L. In the 
ordered phase, the magnetization is extensive [1l oc L, whereas in the disordered phase, the magnetization of clusters 



11 



remain finite /iL oc Cst as it should. At criticality, we measure /ij, « L df where the fractal dimension of critical 
ordered clusters is of order dt ~ 0.85 whereas the exact result is 

d exact = = 0.809... (49) 

This discrepancy will be discussed after Eq. 1571 



F. Critical exponents in the disordered phase h > h c 




-13 -11 



-5 -3 
In (h-h ) 



-13 -11 



-5 -3 
In (h-h ) 



FIG. 4: (d — 1 ) Critical exponents in the disordered phase h > h c : (a) Divergence of the typical correlation length £ tyP of 
Eq. [50] : we measure vt yv ~ 1 (b) Essential singularity of the typical hnite renormalized transverse held /i^ p of Eq. [6] : we 
measure k ~ 1. 



In the disordered phase, the exponential decay of the typical renormalized coupling J^ yp = e ln Jl defines the typical 
correlation length £ typ (EqO 



lnJl yp = lnJ L ~ —f- 
As shown on Fig. 0] (a), the divergence at criticality 

is governed by the exponent 

v typ — 1 



(50) 



(51) 



(52) 



in agreement with the exact solution [3] . 

As shown on Fig. H](b), the finite typical renormalized transverse field h 1 ^ displays the essential singularity (Eq. 
EJ with 



k ~ 1 



(53) 



in agreement with the exact solution j4J . Both values of Eq. [52] and Eq. [53] correspond to the finite-size correlation 
exponent (see Eqs [7] and HP]) 



v FS ~ 2 



(54) 



in agreement with the exact solution 
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FIG. 5: (d = 1 ) Singularity of the intensive magnetization m in the ordered phase (Eq. I14|l : we measure an exponent of 
order f3 ~ 0.33 instead of the exact value of Eq. [57] (see explanations in the text). 

G. Critical exponents in the ordered phase h < h c 

In the ordered phase, the exponential decay of the typical renormalized transverse field defines the length ^ (Eq. 

lnh^ p = Wh^ ~ ~ (55) 

where £/, diverges with the exponent Vh — 1 in agreement with the exact solution 4]. Similarly, the finite typical 
renormalized coupling J^ p displays an essential singularity with the same exponent k ~ 1 as in Eq. 1531 in agreement 
with the exact solution [4] . On Fig. [5] we show our data concerning the intensive magnetization of Eq. [14] 

fn = oc (h c — hf (56) 
We measure an exponent of order (3 ~ 0.33 instead of the exact result 

pexact = ^_ ^5 = 38195 (57) 

We believe that this discrepancy, related to the previous discrepancy found at criticality (see Eq |4"5]) comes from 
the rare-event nature of the magnetization in the critical region : indeed, the critical magnetization represents some 
'persistence exponent' for the strong disorder renormalization flow [25|, because it is related to the probability for 
a given spin to remain undecimated during the RG flow [4]. Within our framework where some spins are declared 
'undecimable' up to some given stage of the real space RG, this persistence probability is changed so that the 
corresponding exponents for the magnetization are not captured exactly by our fixed-cell-size procedure. Nevertheless, 
the approximated exponents obtained are expected to become better and better as the rescaling factor b grows, and 
to converge to the exact values as b — > +oo 

H. Discussion 

In summary, except for the magnetic exponent f3 and related exponents like df that are not reproduced exactly 
for finite b (for reasons explained just above), the fixed cell-size procedure with the rescaling factor b = 2 is able to 
capture correctly the other critical exponents, in particular the activated exponent tp, the typical correlation exponent 
v t y P , and the essential singularity exponent k = 1. For other rescaling factors 6 = 3,4, .., we expect that these critical 
exponents (i/>, vtyp, ft) will be again exactly reproduced (since b = 2 is the 'worst' approximation with respect to the 
full rules corresponding to b — > +oo), whereas the magnetic exponents (3(b) will be better approximations of the true 
exponent f3 exa ct = P(b —> +oo) than (3(b = 2), but will never be exact for any finite b. We also believe that these 
conclusions can be extended to more complicated observables like spin-spin correlations as follows : typical observables 
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(like typical correlations that involve the exponent ip at criticality) should be reproduced with the correct critical 
exponents, whereas averaged observables governed by rare persistent events (like averaged correlations that involve 
the magnetic exponent 0) cannot be correctly reproduced for any finite b. 

After having checked that typical exponents can be correctly reproduced in d = 1 with b = 2, we may now apply 
this fixed cell-size procedure to fractal lattices of dimension d > 1. 



IV. NUMERICAL RESULTS FOR THE SIERPINSKI GASKET 



The Sierpinski gasket is one of the simplest hierarchical lattice made of nested triangles. It has for fractal dimension 

1 3 

D = — = 1.5849525... (58) 
There is no underlying ferromagnetic phase (i.e. T c = 0), because the so-called 'ramification' remains finite (26[. 



A. Principle of the fixed cell-size RG procedure 



A (r A , A(I A ) 




(r B ,A[i B ) (r c ,A(i c ) 

FIG. 6: To construct an open triangle of generation (n + 1) with its renormalized variables 

(Jab, J ac, Jbc,ta, A{ia,tb, Afj,B,rc, A/j,c), the are two steps : 

Step 1 : one draw three independent open triangles of generation n, and one connects them via intermediate points 
Step 2 : one applies the usual strong disorder RG rules to the internal structure. 



Since the Sierpinski gasket can be constructed recursively from triangles, it is clear that the fixed 
cell-size RG procedure based on Strong Disorder RG rules will be based on the joint probability 
Pu(Jab, J ac, Jbc, t a, A/iB,rC) A/ip) of 'open triangles', that will replace the joint probability 

Pti{ Jab, t a, A^; tb, A^b) of 'open bonds' described in section lnB- II : the variables for one open triangle ABC are the 
three ferromagnetic couplings (Jab, J ac-, Jbc) associated to the three bonds of the triangle, the three magnetization- 
excesses (AfiA, A/zs, A/ic) associated to the three vertices, and the three multiplicative factors (ta, tb, tc) for the 
transverse fields of the three vertices. 

It is clear that with three open triangles of generation n, we can build the structure of generation (n + 1) by 
following the same principles as in section III B 21 We may then apply Strong Disorder RG rules to this structure to 
obtain an open triangle of generation (n + 1) with its renormalized variables, by following the same principles as in 
section Til B 31 (the only novelty is that in the ordered phase, one of the three renormalized ferromagnetic coupling of 
the triangle may vanish). 

As in the previous section, we have used the flat initial distribution of Eq. 0D] for the random couplings and a 
uniform transverse field h. The results presented in this Section have been obtained with a pool of size M = 7.10 6 , 
leading to a transition point of order 



1.3527037 < h p ° ol (M = 7.10 6 ) < 1.3527038 



(59) 
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using < n < 100 generations corresponding to lengths 



1 < L = 2™ < 2 100 



(60) 



B. RG flow of the typical renormalized coupling J l ^ v 




70 



In A, 




FIG. 7: (Sierpinski ) (a) RG flow of the logarithm of the typical renormalized coupling (In J^ ap ) in a log-log plot : the slope is 
in the ordered phase (e.g. h = 1), the slope is 1 in the disordered phase (e.g. h — 1.7), and the slope is ip ~ 0.58 at criticality 
(before the bifurcation of the curves corresponding to h — 1.3527037 and h = 1.3527038). (b) RG flow of the width Ax, of the 
distribution of the logarithm of the renormalized couplings in a log-log plot : the slope is in the ordered phase (e.g. h = 1), 
the slope is lo ~ 1 in the disordered phase (e.g. h — 1.7), and the slope is ip ~ 0.58 at criticality (before the bifurcation of the 
curves corresponding to h — 1.3527037 and h — 1.3527038 ). 



On Fig. [7] (a) , we show the RG flow of the typical coupling J 



typ 
L 



; ln J L ag a f uncT j on f Jj 



lnJ^ p 


h<h c 




oc Cst 




L- 




lnjf p 


h=h c 




oc —ifi with V ^ 0.58 




L- 




ln4 w 


h>hc 


L- 


oc — L 



(61) 



C. RG flow of the width of the logarithms of the renormalized couplings 

On Fig. [7] (b), we show the RG flow of the width of the distribution of the logarithms of the couplings 



Ai|ft, <hc oc Cst 

A L \ h=hc oc with V^0.58 

L— > + oo 
L— >+oo 



(62) 



Note that the value oj = 1 is an 'anomalously' large exponent for the droplet exponent of the Directed Polymer if one 
compares with hypercubic lattices where to < 1 (Eq. [5J. We believe that this anomaly comes from the hierarchical 
structure of the Sierpinski gasket where the polymer is forced to visit certain points at each scale. 
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FIG. 8: (Sierpinski ) (a) RG flow of the logarithm of the typical renormalized transverse field (ln.h i ^ p ) in a log-log plot : the 
slope is D = rM; = 1.58.. in the ordered phase (e.g. h = 1), the slope is in the disordered phase (e.g. h = 1.7), and the slope 
is ip ~ 0.58 at criticality (before the bifurcation of the curves corresponding to h = 1.3527037 and h = 1.3527038 ). (b) RG 
flow of the renormalized magnetization fxt in a log-log plot : the slope is D = J^-l = 1.58.. in the ordered phase (e.g. h = 1), 
the slope is in the disordered phase (e.g. h — 1.7), and the slope is df ~ 1.08 at criticality (before the bifurcation of the 
curves corresponding to h = 1.3527037 and h = 1.3527038). 

D. RG flow of the typical renormalized transverse field h£ p 



On Fig. |8](a), we show the RG flow of the typical renormalized transverse field h* p as a function of L 

lnh^ p \ h<hc oc -L D 

L— >+oo 



lnh^ p \ h=hc 



oc -L 4 ' with ip ^ 0.58 

L— »+oo 



\nh^ p \ h>ha oc Cst 

L—>+oo 



(63) 



E. RG flow of the renormalized magnetization hl 

On Fig. [8] (b), we show the RG flow of the renormalized magnetization /xj, of surviving clusters as a function of L 

V>L\h<h oc L D 

L— y+oo 

VL\h=h c oc L df with d f ~ 1.08 

At criticality, the exponent a; of the intensive magnetization of Eq. [5] is thus of order 

x = D-d f ~ 0.5 (65) 

F. Critical exponents in the disordered phase h > h c 



In the disordered phase, the exponential decay of the typical renormalized coupling J l ^ v = e ln Jl defines the typical 
correlation length £ typ (Eq. [9]) As shown on Fig. [9] (a) , the divergence near criticality (Eq. [51]) is governed by the 
exponent 

u typ ^ 0.63 (66) 

As shown on Fig. [9] (b), the finite typical renormalized transverse field displays the essential singularity of Eq. 
[5] with the value 

k ~ 0.82 (67) 
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FIG. 9: (Sierpinski ) (a) Divergences of the typical correlation length £t yp of the disordered phase and of the correlation length 
£,h of the ordered phase : we measure vt yp — 0.63 and v\ ~ 0.91 (b) Essential singularities of the typical finite renormalized 
transverse field in the disordered phase and of the typical finite renormalized coupling J^ p in the ordered phase : we 
measure k ~ 0.82 for both. 



G. Critical exponents in the ordered phase h < h c 




In (h -h) 

c 

FIG. 10: (Sierpinski ) Singularity of the intensive magnetization m in the ordered phase (Eq. I14|) : we measure /3 ~ 0.73 



In the ordered phase, the exponential decay of the typical renormalized transverse field Yi^ 9 — e lntlL defines some 
characteristic length ^ (Eq H5]l . As shown on Fig. [§] (a) , we find that £/j diverges with the exponent 

v h ~ 0.91 (68) 

As shown on Fig. |9](a), the finite typical renormalized coupling J^ p displays an essential singularity with the same 
exponent k ~ 0.82 as in Eq. |67l 

Our data concerning the intensive magnetization of Eq. [14] are shown on Fig [10] : we measure an exponent of order 
^-0.73. 

Our various measures are thus compatible with a finite-size correlation exponent of order 

v FS ^ 1.45 (69) 
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FIG. 11: Hierarchical construction of the diamond lattice of branching ratio k = 2 and rescaling factor b — 8 : we show the 
generation n = with a single bond, the generation n = 1 containing kb = 16 bonds organized in k — 2 branches of b — 8 
bonds in series, and generation n = 2 obtained by iteration. 



V. NUMERICAL STUDY FOR A HIERARCHICAL LATTICE HAVING AN UNDERLYING 

CLASSICAL TRANSITION 

Fractal lattices present a classical ferromagnetic transition only if their ramification is infinite (26j . Exactly renor- 
malizablc lattices with a finite connectivity have a finite ramification and do not present a classical ferromagnetic 
transition, as the Sierpinski gasket discussed in the previous section. Exactly renormalizable lattices with a classical 
ferromagnetic transition are thus hierarchical lattices presenting a growing connectivity, such as the lattice shown on 
Fig. QT]that we study in this section : this lattice is constructed recursively from a single link called here generation 
n = 0. At generation n — 1, this single link has been replaced by k = 2 branches, each branch containing 6 = 8 bonds 
in series. The generation n = 2 is obtained by applying the same transformation to each bond of the generation 
n = 1. At generation n, the length L n between the two extreme sites A and B is L n = b n , whereas the total number 
of bonds grows as N n = (kb) n — so that the fractal dimension reads 

ln(kb) 4 , s 

This type of hierarchical lattices has been much studied in relation with Migdal-Kadanoff block renormalizations 
[27| that can be considered in two ways, either as approximate real space renormalization procedures on hypercubic 
lattices, or as exact renormalization procedures on certain hierarchical lattices [H,|2{|. Various classical disordered 
models have been studied on these hiearchical lattic e, in p articular the diluted Ising model 130 ] , the ferro mag netic 
random Potts model [H , spin-glasses @, S3, and the directed polymer model 0, [H EJSf . 

As in the previous sections, we have used the flat initial distribution of Eq. 001 for the random couplings and a 
uniform transverse field h. The results presented in this Section have been obtained with a pool of size M = 10 7 . We 
find that the corresponding transition point satisfies 

0.556626170 < h p c ool (M = 10 7 ) <h = 0.556626171 (71) 

using < n < 42 generations corresponding to lengths 

1 < L = 8™ < 8 42 ( 72) 

A. RG flow of the typical renormalized coupling J l ^ v 



On Fig. [12j we show our data concerning the RG flow of the typical coupling J^ p = e ln Jl at a function of L 
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(73) 
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FIG. 12: (Diamond k — 2, b — 8) (a) RG flow of the logarithm of the typical renormalized coupling (In J L VP ) in a log-log plot 
: the slope is 1 in the disordered phase (e.g. h — 0.9), and the slope is ip ~ 0.49 at criticality (before the bifurcation of the 
curves corresponding to h — 0.556626170 and h = 0.556626171 ). (b) RG flow in the ordered phase : the asymptotic growth 
corresponds to the slope D a = D — 1 = 1/3 (see Eq. [73}. 

The novelty with respect to the previous cases is thus the asymptotic growth J L VP oc L Ds in the ordered phase, as 
expected when there exists an underlying classical transition (see the discussion around Eq. [18]). 

B. RG flow of the width Al of the logarithms of the renormalized couplings 
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FIG. 13: (Diamond k — 2, b — 8) (a) RG flow of the width Al of the distribution of the logarithm of the renormalized 
couplings in a log-log plot : the slope is uj ~ 0.41 in the disordered phase (e.g. h — 0.9), and the slope is ip ~ 0.49 at criticality 
(before the bifurcation of the curves corresponding to h — 0.556626170 and h = 0.556626171 (b) In the ordered phase, the 
width Al decays asymptotically as L LJ s- D s with uj s — D s ~ —0.29 (see Eq. I74p 

On Fig. [T31 we show our data concerning the RG flow of the width of the distribution of the logarithms of the 
couplings 



&h\h<h c 
^■h\h=h a 
A L \h>h c 



CX L 
L— y+oo 



(X 
L— >+oo 



Ds with lo s — D s 
with V ^ 0.48 
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oc L w with uj ~ 0.41 
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(74) 
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Here the novelty with respect to the previous cases is again the ordered phase with the decay of the width Al oc L Ws_Ds 
as expected when there exists an underlying classical transition (see the discussion around Eq. IT5|) . 

Another important result is the inequality ip > u> implying a singular diverging amplitude A for the width in the 
disordered phase (see the discussion around Eq. [T2|) . 
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FIG. 14: (Diamond k — 2, b = 8) (a) RG flow of the logarithm of the typical renormalized transverse field (In h 1 ^) in a 
log-log plot : the slope is D — 4/3 in the ordered phase (e.g. h — 0.2), the slope is D s = 1/3 in the disordered phase (e.g. 
h = 0.9), and the slope is ip ~ 0.49 at criticality (before the bifurcation of the curves corresponding to h — 0.556626170 and 
h = 0.556626171 ). (b) RG flow of the renormalized magnetization /ii in a log- log plot : the slope is D = 4/3 in the ordered 
phase (e.g. h — 0.2), the slope is D s — 1/3 in the disordered phase (e.g. h — 0.9), and the slope is df ~ 0.85 at criticality 
(before the bifurcation of the curves corresponding to h — 0.556626170 and h = 0.556626171). 



C. RG flow of the typical renormalized transverse field h^ p 

On Fig. [TH (a), we show the RG flow of h l ^ p as a function of the RG scale L 

L D 

Lfi with ifj ~ 0.49 

L D ° with D S = D- 1 = 1/3 (75) 

The novelty with respect to the previous cases is the disordered phase that involves the dimension D s as a consequence 
of the growing connectivity of the diamond hierarchical lattice (see the discussion in I V Gp . 

D. RG flow of the renormalized magnetization [il of surviving clusters 

On Fig. [T3] (b), we show the RG flow of the renormalized magnetization (Xl of surviving clusters as a function of 
the RG scale L 

L— > + oo 

fiL\h=h c « L<is with d f ~ 0.85 

L^-\-oo 

fiL\h>h e « L Db with D S =D- 1 = 1/3 (76) 

L— >+oo 

At criticality, the exponent x of the intensive magnetization of Eq. [S] is thus of order 

x = D-d f - 0.48 (77) 



lah t ^ p \h<h e oc - 

L—>-\-oo 

Inh^lh^ ex - 
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\n hf p \ h> h c oc - 
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E. Critical exponents in the disordered phase h > h c 




FIG. 15: (Diamond k = 2, b = 8) Critical exponents (a) Divergences of the typical correlation length £t yP of the disordered 
phase and of the correlation length £h of the ordered phase : we measure v tyv — 0.86 and v\ — 1.08 (b) Singularity of the 
intensive magnetization m in the ordered phase (Eq. I14[) : we measure ~ 0.81. 



In the disordered phase, the exponential decay of the typical renormalized coupling J^ p = e L defines the typical 
correlation length £ typ As shown on Fig. [15] (a), the divergence near criticality is governed by an exponent of order 

Vty V = °- 86 (78) 



F. Critical exponents in the ordered phase h < h c 



In the ordered phase, the exponential decay of the typical renormalized transverse field h£ p — e lnllL defines the 
typical cluster linear length £/, 

kfc^ETir * -(pj D (79) 

L^ + oo \t,hj 

where ^ diverges with the exponent Vh — 1.08 (see Fig [TBI (a) - ) 

On Fig. [TS] (b), we show our data concerning the intensive magnetization : we measure an exponent of order 
/3 ~ 0.81. 

Our various measures are thus compatible with a finite-size correlation exponent of order 

v FS ^ 1.71 (80) 



G. Discussion : effect of the growing connectivity 



In this section, we have considered a hierarchical fractal lattice, where the end points have a growing connectivity 
as k n , when the length grows as L n = b n . This property has for consequence that in the disordered phase, the 
magnetization /i^ and the logarithms of the transverse fields do not remain finite as they do when they represent sites 
of finite connectivity, but scale as the growing connectivity as k n = L^ s where 

„ In k „ , „ 

= hib = ( 81 ) 

By consistency in the disordered phase, only sites should be decimated asymptotically, i.e. the logarithms of the 
renormalized transverse fields scaling as 



ln/ii cx — L 

L— >+oo 



(82) 
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should remain bigger than the logarithm of the renormalized couplings scaling as in Eq. [8] : the RG procedure is thus 
consistent only if Dg < 1. In addition, if one wishes that the spurious behavior of Eq. 155] docs not affect the critical 
point either, one should have the stronger inequality 

D s < r/> (83) 

This condition is satisfied in the case we have studied numerically with D s ~ 0.33 < 0.49 ~ ip (and this is why we 
have chosen the value 6 = 8 instead of smaller values). But it is clear that this condition will not be satisfied by all 
diamond hierarchical lattices with arbitrary values of the parameters (fc, b) of Eq. 1701 



VI. CONCLUSION 



In this paper, we have proposed to include Strong Disorder RG ideas within the more traditional fixed cell-size real 
space RG framework. We have first considered the one-dimensional chain as a test for this fixed cell-size procedure. 
Our conclusion is that all exactly known critical exponents are reproduced correctly, except for the magnetic exponent 
/3, because it is related to more subtle persistence properties of the full RG flow. We have then applied numerically 
this fixed cell-size RG procedure to two types of renormalizable fractal lattices : 

(i) the Sierpinski gasket of fractal dimension D — In 3/ In 2, where there is no underlying classical ferromagnetic 
transition, so that the RG flow in the ordered phase is similar to what happens in d = 1. 

(ii) a hierarchical diamond lattice of fractal dimension D = 4/3, where there is an underlying classical ferromagnetic 
transition, so that the RG flow in the ordered phase is similar to what happens on hypercubic lattices of dimension 
d> 1. 

In both cases, we have found that the transition is governed by an Infinite Disorder Fixed Point. Besides the measure 
of the activated exponent ip at criticality, we have analyzed the RG flow of various observables in the disordered phase 
and in the ordered phase, in order to extract the 'typical' correlation length exponents of these two phases (respectively 
vtyp of Eq. IH1 and v\ x of Eq. II 5|) which arc different from the finite-size correlation exponent vps that governs all 
finite-size properties in the critical region (Eq. [3]). In the disordered phase, we have also measured the fluctuation 
exponent w (Eq. [5]) which is expected to coincide with the droplet exponent of Directed Polymer living on the same 
lattice (20l |. Whereas the two exponents t/> and u> are known to coincide in dimension d = 1, we have found here in 
section IVl that it is possible to have oj < ip (and accordingly a diverging amplitude A(h) in Eq. |5J. Moreover the 
measured value ip ~ 0.49 is close to the values measured for hypercubic lattices in dimensions d = 2, 3, 4 |6l-[l6j|. 

Besides the random transverse field Ising model that we have discussed, we hope that the idea to include Strong 
Disorder RG principles within a fixed cell-size real space renormalization will be useful to study other types of 
disordered models on fractal lattices, like for instance random walks with disorder on various fractal lattices that have 
been studied recently in [49l |. 



Appendix A: Reminder on the full Strong Disorder RG procedure in energy 

In this section, we recall the standard Strong Disorder Renormalization for the Random Transverse Field Ising 
Model of Eq. [U to compare with the fixed cell-size framework introduced in the text. 



1. Reminder on Strong Disorder RG rules on arbitrary lattices 

For the model of Eq. [TJ the Strong Disorder RG rules are formulated on arbitrary lattices as follows [f| 0] : 
(0) Find the maximal value among the transverse fields hi and the ferromagnetic couplings Jjk 

£1 = max [hi, Jjk] (Al) 

i) If Q, = hi, then the site i is decimated and disappears, while all couples (J, k) of neighbors of i are now linked via 
the renormalized ferromagnetic coupling 

J jk - Jjk H ^— (A2) 

ii) If Q — Jij, then the site j is merged with the site i. The new renormalized site i has a reduced renormalized 
transverse field 

hT w = h in with n = (A3) 
•Jij 
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and a bigger magnetic moment 

f iT = tH +H j (A4) 
This renormalized cluster is connected to other sites via the renormalized couplings 

JIT = J >k + Jjk (A5) 

(hi) return to (0). 

2. RG flow of probability distributions 

The above RG rules determine the RG flow for the joint probability distribution at scale 

Q = er r (A6) 
of the renormalized clusters characterized by their magnetic moments /i£, their renormalized fields 

hi = tte- ft = e-( r +*) (A7) 
and the renormalized couplings between them 

= fie-* 1 * = e"( r+c «^ (A8) 

where ft and Qj are positive random variables. In the following, we recall 0, uj the properties of the probability 
distributions i?r(ft) and iV(C*i) (normalized per surviving cluster) and of the density nr of surviving clusters per 
unit volume that defined the characteristic length-scale lr via 

nr ~ 4 (A9) 
r 

3. Disordered phase 

In the disordered phase, renormalized clusters remain 'finite', i.e. asymptotically at large RG scale T, only transverse 
fields are decimated via the rule (i) of section |A~T1 while the rule (ii) does not occur anymore. This means that the 
variable ft of Eq. IA7I remains finite, i.e. the probability distribution i?r(ft) converges towards a finite probability 
distribution i?oo(ft) without any rescaling for the variable ft. The value -Roo(O) at the origin (3 = will then govern 
the decay of the clusters density nr 

d r n r = -Roo(0)n r (A10) 

leading to the exponential decay 

nr«e- rJi »M (All) 



or equivalently to the exponential growth of the length-scale lr of Eq. IA9I representing the typical distance between 
surviving clusters 

Hoa (0) 

lr oc e r ^— (A12) 
The relation between the energy scale f2 = e _r and the length scale lr thus corresponds to the power-law 

n = e~ r cx Zf z (A13) 



with the continuously variable dynamical exponent 
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4. Critical point 

At the critical point, both transverse fields and couplings continue to be decimated at large scale T, and there is no 
finite characteristic scale : as a consequence, the only scale available for the variables /3j of Eq. I A7I and Qj of Eq. IA8I 
is the RG scale r itself. More precisely, the probability distributions Pr(A) an d Pr(Cij) (normalized per surviving 
cluster) follow the scaling form 













Cij 




» r 



RrWi) 



where the rescaled probability distributions 72. and P do not depend on T. 

The values 72(0) and P(0) at the origin will then govern the decay of the clusters density nr 



"r oc rn0 ) +no) ( Al7 ) 



leading to the power-law 



or equivalently to the power-law growth of the length-scale lp of Eq. IA9I 

l r oc r * (A18) 
The relation between the energy scale ft — e~ T and the length scale It thus corresponds to the activated form 

n = e~ r cx e~'r (A19) 

where the activated exponent ip reads 

* = 72(0) +P(0) (A20) 

Finally, the magnetic moment of surviving clusters is expected to grow as a power-law of T, and equivalently as 
some power-law of Zr (Eqs IA18I and IA20|I 

/iocr*« with d f = <S>iP (A21) 
where df represents the fractal dimension of surviving critical clusters. 

5. Ordered phase 

In the ordered phase, renormalized clusters are expected to become 'extensive', i.e. asymptotically at large RG 
scale r, only couplings are decimated via the rule (ii) of section[ATJ while the rule (i) does not occur anymore. For the 
asymptotic behavior of the renormalized couplings between surviving clusters, one has to distinguish two possible 
cases, depending on the existence or non-existence of an underlying classical ferromagnetic phase. 

a. Ordered phase when there is no underlying classical ferromagnetic phase (as in d — 1) 

When there is no underlying classical ferromagnetic phase (as in the one-dimensional chain), the variables Qj of 
Eq. I A8 1 remain finite, i.e. the probability distribution Pr{Cij) (normalized per surviving clusters) converges towards a 
finite probability distribution (£y ) without any rescaling for the variable C,ij ■ The value (0) at the origin £ = 
will then govern the decay of the clusters density nr 

drriT = -Poo (0)nr (A22) 
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leading to the exponential decay 

n r oce- rp ~(°) (A23) 

or equivalently to the exponential growth of the length-scale lr of Eq. IA9I representing the typical linear size of 
surviving clusters 

l r oc e r ^ (A24) 

The relation between the energy scale O = e~ r and the length scale lr thus corresponds to the power-law of Eq. IA13I 
with the continuously variable dynamical exponent 

(A25) 



Poo(0) 



i.e. the ordered phase is thus rather 'symmetric' to the disordered phase (see Eq IA14|) . Note that in dimension d = 1, 
this symmetry is a consequence of the duality between couplings and transverse fields Q. 



b. Ordered phase when there is an underlying classical ferromagnetic phase (as in d > 1) 

When there exists an underlying classical ferromagnetic phase (as on the hypercubic lattices in dimension d > 1), 
the ordered phase is completely different from the properties just described, because an infinite percolation cluster 
appears at some finite RG scale T perco : we refer to Ref Q for more details. 
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